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Abstract 

We study 1-dimensional coupled map lattices consisting of diffusively coupled Cheby- 
shev maps of A^-th order. For small coupling constants a we determine the invariant 
1-point and 2-point densities of these nonhyperbolic systems in a perturbative way. For 
arbitrarily small couplings a > the densities exhibit a selfsimilar cascade of patterns, 
which we analyse in detail. We prove that there are log-periodic oscillations of the den- 
sity both in phase space as well as in parameter space. We show that expectations of 
arbitrary observables scale with y/a in the low-coupling limit, contrasting the case of hy- 
perbolic maps where one has scaling with a. Moreover we prove that there are log-periodic 
oscillations of period logA^^ modulating the y^-dependence of the expectation value of 
any given observable. 
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1 Introduction 



Coupled map lattices (CMLs) as introduced by Kaneko and Kapral [H [2] are interesting exam- 
ples of spatially extended dynamical systems that exhibit a rich structure of complex dynamical 
phenomena [21 HI [5l [H [7j . Both space and time are discrete for these types of dynamical systems, 
but the local field variable at each lattice site takes on continuous values. CMLs have a variety 
of applications: They have been studied e.g. as models for hydrodynamical flows, turbulence, 
chemical reactions, biological systems, and quantum field theories (see e.g. reviews in [31 II])- 
Of particular interest are CMLs that exhibit spatio-temporal chaotic behaviour. The analysis 
of these types of CMLs is often restricted to numerical investigations and only a few analytical 
results are known. A notable exception is the case of CMLs consisting locally of hyperbolic 
maps for small couphng a. Here a variety of analytical results is known [8], [HI [101 [H] that 
guarantee the existence of a smooth invariant density and ergodic behaviour. 

The situation is much more complicated for nonhyperbolic local maps which correspond to 
the generic case of physical interest [121 [131 [IH [HI [HI [T71 [18]. In the one-dimensional case, 
non-hyperbolicity simply means that the absolute value of the slope of the local map is allowed 
to be smaller than 1 in some regions. For chaotic CMLs consisting of nonhyperbolic maps 
much less is known analytically, and standard techniques from ergodic theory are usually not 
applicable. From the applications point of view, the nonhyperbolic case is certainly the most 
interesting one. For example, it has been shown that nonhyperbolic CMLs exhibiting spatio- 
temporal chaos naturally arise from stochastically quantised scalar field theories in a suitable 
limit where the strength of the underlying self- interacting potentials goes to infinity |H [TU] . 
CMLs obtained in this way can serve as useful models for vacuum fluctuations and dark energy 
in the universe [20]. Other applications include chemical kinetics as described by discretised 
reaction diffusion dynamics [3]. 

In this paper, we study weakly coupled A^-th order Chebyshev maps [HI [191 [201 [211 [221 [23] • 
We consider CMLs with diffusive coupling, i.e. there is a coupling between neighboured maps 
that resembles a discrete version of the Laplacian on the lattice. This form of the coupling is 
the most important one for the quantum field theoretical applications [1]. In the uncoupled 
case it is well known that the Chebyshev maps are conjugated to a Bernoulli shift of symbols 
(see, e.g. [21] )■ In the coupled case, this conjugacy is destroyed and the conventional treatments 
for hyperbolic maps does not apply, since the Chebyshev maps have A^ — 1 critical points where 
the slope vanishes, thus corresponding to a nonhyperbolic situation. 

We will derive explicit formulas for the invariant densities describing the ergodic behaviour 
of weakly coupled Chebyshev systems, by applying a perturbative treatment. Our analytical 
techniques will yield explicit perturbative expressions for the invariant 1-point and 2-point 
density for small couplings a. We will prove a variety of interesting results concerning the 
scaling behaviour of these nonhyperbolic systems. Our main results are that the 1-point density 
exhibits log-periodic oscillations of period log A^^ near the edges of the interval. Expectations of 
typical observables scale with y/a (rather than with a as for hyperbolic coupled maps). We also 
show that there are log-periodic modulations of arbitrary expectation values as a function of 
the coupling parameter a. Although our results are rigorously derived for the explicit example 
of Chebyshev maps only, we expect that similar results hold for other nonhyperbolic systems 
as well, since our techniques are quite general and can be applied to various other systems in 
a similar way. 
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This paper is organized as follows: In section 2 we introduce the relevant class of coupled 
map lattices and provide some numerical evidence for the observed scaling phenomena, choosing 
as a main example the Ulam map (the N = 2 Chebyshev map). Our perturbative theory of 
the invariant density of the coupled system is presented in section 3. Based on these results we 
explain details of the selfsimilar patterns in section 4. In section 5 we extend our results to the 
third order Chebyshev map, where, due to an additional discrete symmetry, the behaviour is 
somewhat different as compared to iV = 2. Finally, in section 6 we collect our results to finally 
prove that there are log-periodic oscillations modulating the scaling behaviour of expectation 
values of arbitrary observables if the couphng parameter a is changed. 

2 Observed scaling phenomena 

We consider one-dimensional coupled map lattices of the form 

0^+1 = (1 - a)T^m + + TU<PI-')) . (1) 

Here T/v(0) denotes the A^-th order Chebyshev polynomial and 0^ is a local iterate at lattice 
site i at time n. Our main interest is concentrated on the cases N = 2 where T2(0) = 20^ — 1 
and = 3 where T^{(f)) = 40^ — 30. The sign s = ±1 distinguishes between the 'diffusive 
coupling' [s = +1) and the ' ant i- diffusive coupling' {s = —1). The index b finally can take 
the values 6 = 0, 1 where 6 = 1 corresponds to forward coupling (T^(0) := T/v(0)) whereas 
6 = stands for backward coupling (T^(0) := 0). The discrete chaotic noise field variables 0^ 
take on continuous values on the interval [—1, 1]. The initial values 0q are randomly chosen in 
this interval with a smooth density, usually chosen to be the invariant density of the uncoupled 
system. 

2.1 Scaling in parameter space 

The variable a in Eq. ([1]) is a coupling parameter which takes values in the interval [0, 1]. It 
determines the strength of the Laplacian coupling in Eq. ([T]) (in the particle physics applications 
described in |1] can be regarded as a kind of metric for a one-dimensional 'chaotic string'). 
For a = we end up with the chaotic behaviour 0^+i = T/v(0^) as generated by uncoupled 
Chebyshev maps. If a small coupling a > is switched on, one observes the formation of 
nontrivial density patterns with scaling behaviour with respect to y/a. Consider an arbitrary 
test function (observable) V{(f)) of the local iterates, and let (. . .)a denote the expectation values 
formed by long-term iteration and averaging over all lattice sites. For small a one observes 
numerically 

(V(0)),-(y(0))o = /(^)(lna)aV2 (2) 

where /*^^)(lna) is a periodic function of In a with period InA^^. The function f^^^ = f^'^^[V] 
is a functional of V. A double logarithmic plot of |(V^(0))a — (^(0))o| versus a hence shows 
a straight line that is modulated by oscillating behaviour. Examples are shown in Fig. 1. We 
will rigorously prove the observed log-periodic scaling behaviour with respect to a in section 6 
of this paper. 

For small a there is no difference between expectations calculated with the forward or 
backward coupling form of ([1]), the value of 6 matters only for large values of the coupling. 
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Figure 1: Scaling in parameter space for the test functions V{(f)) = V2(0) := — (A^ = 2) 
and ^(0) = V3(0) := — 0"^ (A^ = 3). To calculate the expectations, the coupled maps were 
iterated 10^ times for random initial conditions. The lattice size was 10^. 



2.2 Scaling in the phase space 

There is also scaling behaviour in the phase space of the coupled system. First, let us consider 
the uncoupled case. The invariant (one-point propability) density corresponding to a = is 
known to be 



(see e.g. for a derivation). This density is universal in the sense that it does not depend on 
the index N of the Chebyshev polynomial for N >2. 

The 1-point invariant density at each lattice site will of course change if we consider values 
a 7^ 0. It then depends not only on a but also on the index A^. In other words, the A^- 
degeneracy is broken for finite coupling. In principle, the invariant densities of coupled map 
lattices can be understood by finding fixed points of the Perron-Frobenius operator of these 
very high-dimensional dynamical systems fiUi [25| [26| [271 [28] . The 1-point density can then 
be obtained as a marginal distribution, by integrating out all lattice degrees of freedom up to 
one. In practice, this program is rather complicated for our type of systems and we will apply 
a more pedestrian but feasable method in the following sections. 

If the coupling a is slightly increased starting from a = 0, the invariant 1-point density Pa(0) 
starts to deviate from po{(f)). There are particularly strong deviations at the boundaries = ±1 
of the interval. While for odd A^ the evolution equation ([1]) is invariant under the replacement 
(j) —>■ —(f) and hence the invariant density is symmetric, this is not the case for even values of 
A^. In both cases we are interested in the region close to = ±1 and our goal is to understand 
the deviation of Pa{4>) from po{(p). 
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Figure 2: Energy difference Ae(0) for the 2B-string at a = 0.00007625 



For our purpose it is useful to define an 'effective energy' Ca as follows: 

1 



ea(0) 



1 
2 



27r2pa(0)^ 



(4) 



This function is a useful tool if one wants to apply generalized statistical mechanics techniques 
to coupled Chebyshev maps, see [1] for more details. For a = the above effective energy is 
just given by the 'kinetic energy' expression |(/)^ of a free particle, and the invariant density 
coincides with that of a generalized canonical ensemble in Tsallis' formulation of nonextensive 
statistical mechanics ^29j, the entropic index q being given by g = 3. For a 7^ the effective 
energy becomes more complicated and it includes the effects of interactions with neighboured 
lattice sites. The quantity of interest is the energy difference Ae(0) = ea(0) — eo(0). 

In [1] a general classification scheme of the relevant Chebyshev CMLs was introduced. To 
be specific, let us here consider the coupled Chebyshev map system with = 2, s = +1 and 
b = (called 2B-string in [4j). For finite a we obtain an (asymmetric) energy difference Ae(0) 
which when approaching the boundary = +1 shows selfsimilar behaviour, i.e. a repetition of 
nearly the same pattern if a suitable rescaling is done (see Fig. [2]). 

The scaling behaviour close to = +1 can be analysed in more detail by looking at a 



logarithmic plot of Ae(0) versus 



log4(l — 0), as shown in Fig. [31 In order to deal with a 



logarithmic histogram plot the data are adjusted in the usual way by the factor 1/(1 — 0) ln4, 
as obtained from transformation of random variables. We see that the pattern in this plot is 
repeated with a period A0* = 1, though there is a sharp breakdown of periodicity at the value 
0* = a* = — log4(4a). We will prove this behaviour later. 



5 




012345678 

<P* = -log4( 1 

Figure 3: Energy difference Ae(0) for tlie 2B-string for a = 0.00007625 (logaritfimic plot) 



2.3 Connection between parameter and phase space scaling 

Let us now look at liow tlie energy difference cfianges witli tlie coupling parameter a. If the 
coupling is decreased by a factor of 4 one observes that the pattern in the logarithmic plot is 
moved to the right by A0* = 1 while a new pattern element is created to the left. The picture 
is like a wave generator located at = that was started at some 'time' in the past and is 
sending out waves with a speed A0* = 1. In this analogy 'time' corresponds to the negative 
logarithm of the coupling parameter. Since in Fig. [3] the pattern amplitude becomes smaller 
and smaller towards higher values of 0* we can try to compensate this effect by multiplying 
Ae(0) by an appropriate simple function of 0*. It turns out that the multiplication by po(0) 
compensates the decrease of the pattern for a fixed value of the couphng a. However, the 
pattern as a whole turns out to decrease by a factor of 2 if we divide the coupling by a factor 
of 4. Therefore, as a last adjustment, we multiply with the inverse square root l/i/a of the 
coupling. In a sequence of steps Aa* = 1 for the 'time variable' a* = — log4(4a) we show 
in Fig. m the progressive motion of the adjusted pattern, indicating the breakdown point at 
0* = a*. 

We can compensate the shift of the pattern to the right by subtracting a* = — log4(4a) 
from 0* = — log4(l — 0) obtaining the new variable x* = 0* — a*. Note that x* is negative 
beyond the breakdown region. In writing subsequent curves into the same diagram we obtain 
an entire set of curves shown in Fig. [5l While the curves diverge for large absolute values of x*, 
the surprising fact is that they exactly coincide close to the breakdown at x* = 0. This means 
that there is invariance under a suitable renormalization transformation: a~^/^po(0)Ae(0) is a 
function of x* = — log4(a;) only if the absolute value of x* is small. The shape of the first peak 
(i.e. for X* G [—1, 0]) seems to be very simple, while the following peaks increase in complexity. 

As a last diagram illustrating this behaviour, in Fig. [6] we show the development of the 
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Figure 4: Adjusted energy difference a~^/^po (0)^^(0) for tlie 2B-string at a = 0.00007625, a = 
0.0000190625, and a = 0.000004765625 (from top to bottom). Indicated is tlie breakdown point 
of tfie pattern at 0* = a* = — log4(4a) (vertical line) as well as the pattern amplitude multiplied 
by a factor of 0.001 (horizontal hue) in order to have a suitable scale for the breakdown region. 
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Figure 5: Adjusted energy difference a"^^Vo(0)Ae(0) for tlie 2B-string at a = 0.00007625, 
a = 0.0000190625, and a = 0.000004765625 (from top to bottom) as a function of x* ^ (p* - a*. 
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Figure 6: Comparison of tlie first tliree peaks, as obtained for a fixed interval of 0* witli 
a = 0.00007625, a = 0.0000190625, and a = 0.000004765625. The indicated logarithms in the 
diagram are related to the value a — 0.00007625. 
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Figure 7: Invariant density Pa{4>) close to = +1 for the 2B-string for a=0. 00007625, as 
compared to the invariant density po(0) for vanishing couphng (dotted curve). 



density pattern for a given fixed interval of 0* for the already used sequence of values of a. 
Note the slight shift of the maxima to the left and the formation of plateaus. Various details 
of this figure will become clearer by our theoretical treatment in section 4. 

3 Perturbative calculation of the invariant density 

Let us now proceed to an analytic (perturbative) calculation of the invariant 1-point density of 
the 2B-dynamics 

'/'n+l = (1 - «)T2(0^) + ^(0^+^ + (5) 

We are particularly interested in the behaviour near the boundaries = ±1. The numerically 
determined density near = +1 is shown in Fig. [71 

As a first approximation, we essentially regard the system ([5]) as a perturbed one-dimensional 
system, the perturbation of strength a given by the arithmetic mean of nearest neighbours 
{(f)l^^ + 4>l^^)/2 =: (0+ + 0-)/2 =: 0. As a first approximation we may also regard 0+ and 0_ 
to be nearly independent for small a, provided we are not too close to the edges of the interval. 
This assumption is necessary as a starting point for our perturbative approach, corrections 
to this assumption are expected to appear at higher orders of the perturbation series. The 
probability density paa of is obtained by the convolution 

Paai4>) = 2 / Pa(20 - 0__)pa(0-)rf0-, (6) 

the factor 2 arising from the fact that d(j)^d(j)^ = 2(i0(i0_. If we replace the densities pa 
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Figure 8: Invariant densities Poo(0) (solid line) and po(0) (dashed line) 



occurring in Eq. (EI) by po, "we can calculate the convolution analytically to obtain 

2#- '^r.f.r. — 72, 



POO(0) 



l_(2^-^_)Vl- 



4a- 



(7) 



where i^'(a;) is the complete elliptic integral of the first kind and 6{x) is the step function. The 
functions po(0) and poo(0) are shown in Fig. [HI If we calculate to leading order, we can replace 
Pa(0) and Paa(0) by po(0) and poo(0), respectively. For our perturbative approach we always 
assume that the difference between the perturbed and unperturbed density is small, except at 
the edges of the interval. 



3.1 The zeroth iterate of the convolution formula 

Having the density Paa{<P) at hand, we can now calculate the probability that the value (j) 
in an interval of length d(j) and the value in an interval of length d(j) give rise to a value 
(/)' = [1 — a)T2{(j)) + a(j) in an interval of length d(j)', 

Pa(0') = I Pai<P)Paa ((0' - (1 - a)U<P)) / a) dcj). (8) 

Note that there is an implicit constraint in this equation that restricts the integration range 
of (p, given by the condition |0| < 1. So far the equation is quite generally valid. Now let us 
assume that we are close to the left edge, by writing (p' = ay — 1. In this case the integration 
range for (p is restricted by the demand |0| < 1 to a small symmetric interval around 0. We 
may thus write Pa(0) ~ Pa(0) to obtain the leading order approximation of the 1-point density 
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Figure 9: Invariant density Pa{4>) close to = — 1 for the 2B-string with a = 0.00007625, as 
compared to the invariant density po(0) for vanishing couphng (dashed curve) 
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Figure 10: Zeroth iterate pf\(t)) close to = —1, calculated according to Eq. ( fT3l) 
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Paiay — 1) close to (f> = — 1, 



p(^\ay-l) = P-M [\,,{{ay-l-{l-a)T,mici)dct> 



a J -I 

2pa(0) 



a Jo 



Paa {{ay - 1 - (1 - a)r2(0))/a) dcj) (9) 



which we call 'zeroth iterate' (of our convolution scheme). Finally, we can change to an inte- 
gration over 0, 



i(ay-l-(l-a)(202-l)) ^ 0= 



—a ( 



2 J2a(l - a) 



(10) 



to obtain 



pf{ay - 1) = PdM^. (11) 



2a 1 



Note that the integration range in this formula is given by the domain of the square root 
appearing in the denominator and the demand |0| < 1. In the following, we no longer write 
down the integration limits explicitly. Using Eq. ([7]) for the density of 0, we can substitute 
0+ = 20 — 0_ to end up with the symmetric result 

(0)/.. ,^ _ Pa(0) f pa{(t)+)d(t)+pa{(p-)d(t)- 



v/2a(l - a) J 7y-l + (0+ + 0_)/2 

The two integration limits in this formula are given by the condition that the argument of the 
square root should always be positive and that 0+, 0_ G [—1,1]. Note that in calculating the 
leading order approximation of the density close to = —1 there is no need to keep the indices 
a on the right hand side of Eq. f|T2l) . Therefore, we can use the explicit representations in 
Eqs. ([3]) and ([7]) to calculate the integrals. Another possibility is to calculate the integrals as 
time averages. We may just iterate with uncoupled Chebyshev maps T2 to obtain a weighting 
with pq. From this we obtain 

,(0)(a, - 1) = , ' Y. ^'f-^ + '^"^+^"-'/^) (13) 

7rJ2a(l - a) J2 _ 1 + (0^^ + 0„_)/2 



where J is the total number of iterations. The analytic perturbative result calculated in this 
way is shown in Fig. [TUl It coincides quite perfectly with the numerically obtained density 
(obtained by direct iteration of the CML) as shown in Fig. O 
The main conclusions that we can draw up to now are: 

• The behaviour at the boundary close to = —1 is reproduced by the zeroth iterate of 
our convolution scheme with high accuracy. To leading order in a we have Pa{ciy ~ 1) = 

P^^Kay-i). 



12 



Figure 11: Two-point density function PaaiO'V ~ 1;0) in the range 
< 2/ < 3 (horizontally) und — 1 < < +1 (vertically). 



• For arbitrarily small a the density Pa(0) does not approach infinity for — » — 1 (as in the 
uncoupled case) but falls off to after a maximum at = a — 1. 

• If we parametrize the region close to = —1 by = ay — 1, the only dependence on the 
coupling a that is left for the density is Pa{ciy — 1) ~ a^^/^. 

Actually, the last property is already a property of Po(0)) 



However, in combination with the appearence of a last maximum at = a — 1, we obtain a 
scaling relation between a and 0, which induces a variety of interesting results (see section 6). 

3.2 Two-point function and backstep transformation 

The next iteration step of our scheme leads from the region just considered to the main scaling 
region close to = +1 which we parametrize by = 1 — ax. In order to accomplish this step 
we could think of again calculating the convolution according to Eq. ([8]), only that in this case 
we use 0' = 1 — ax. However, for values close to the boundaries neighboured field variables 
exhibit strong correlations, they cannot be regarded as nearly independent anymore. The 2- 
point density function describing the joint probability of neighboured lattice sites exhibits a 
rich structure, as shown in Fig. [TTJ This structure at the border of the interval represents 
the main difficulty to analyse the invariant densities of nonhyperbolic CMLs. Our basic idea 
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(14) 





13 



in the following is to proceed by one backward iteration to the region close to = where 
approximate independence holds, and then transform the densities accordingly. 

In order to construct the two-point function analytically, we make use of the fact that the 
integral of this function paa{ay — 1, 0) over cj) has to reproduce Pa{o,y — 1) since this is just the 
marginal distribution. In searching for a transformation which keeps the integral in Eq. f|T2|) 
invariant but changes the integrand we have to think about an appropriate substitution. The 
one which fits in for our purposes is the transformation = T2^{(p'±) which we may call the 
backstep (or pre-image) transformation. The backstep transformation keeps the unperturbed 
density invariant, 

Po(0±)c?0± = Po{TNi^±))dT^{<P±) = iVpo(0±)#±. (15) 
Under this transformation we obtain quite generally 



where the two-fold sum runs over all N solutions of the inverse map T^^{(j)). For = 2, we 
sum over the two branches of each of the two square roots. In substituting finally 0^ = 20 — 0'_ 
we obtain pf\ay — 1) = / pf}{ay — 1, 0)(i0 where 

(o).„. 1 v-v- /■ 2po(20-0'_)po(0-)rf0'_ 



7r^2a(l-a)22 ± ± J _ i + (±^0 + (i _ 0'_)/2 ± ^(1 + 0'_)/2)/2 

(17) 

This perturbative result reproduces the numerical 2-point distribution observed in Fig. [TTl 
Finally, as a one-point distribution the function p^^}{ay — 1,0) has to be normalised. Actually, 
we obtain this normalisation by no effort because obviously the ratio p^aai^V ~ 1; 0) = P^ai'^V ~ 
l^<p) / pf\ay — 1) (describing a conditional probability) is normalised. 

3.3 The first iteration step: from = — 1 to = +1 

The normalisation is skipped again if we use the conditional probability p^adi'^y ~ 0) close to 
= — 1 for the calculation of the next iterate, 



P^:\4>') = I p^^}{ay-l,{<P' -{l-a)T,{ay-l))la)pf{ay-l)dy 
Pad {ay - 1, (0' - (1 - a)T2{ay - I)) /a) dy. 



For 0' = 1 — ax the second argument reads = 1 + 4y — x + 0{a). Because of this, the range 
of integration for y does not exceed (1 — x)/4 by more than one fourth. Instead of integrating 
over y, we can integrate over which enables us to reverse the backstep transformation, 

pW(l -ax) = ^J Paa (^{X - 1 + 

1 j^j^ f 2po(20-0'_)po(0'„)rf0rf0'_ 



47r J2a(l - a)2^ ^ J{x-1 + 0)/4 - 1 + {T2-\2<j) - 0'_) + T2\<P'_))/2 
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Figure 12: 1-point density pa(l — O'X) at the right edge of the interval (sohd curve) for the 2B- 
string with a = 0.00006725 as a function of x, compared to the first iterate p'^\l — ax) (dashed 
curve), the second iterate p^aK^ ~ '^^) (dotted curve), and the sum of both (dashed-dotted 
curve) 



1 j^T. Po{^'+)d4>'^po{4>'_)d^'^ 



An^2ail - a)2^ J ^{x~l + {<P'+ + 0'_)/2)/4 - 1 + iT,-\<P'+) + T^\<P'.))/2 

1 r poi(l)+)d(f)+po{(f)_)d(t)^ 



47r^2a(l - a) J ^(x - 1 + (T2(0+) + T2(</._))/2)/4 - 1 + (0+ + 0_)/2 

1 f Po(0+)c?0+Po(0-)#- 



47ry'2a(l - a) ^x/A + (0+ + 0_ - 2)/2 + (T2(0+) + T2(0_) - 2)/8 ' 



(19) 



The sums over the branches disappears again. For the last line we have used an equivalent 
representation which will later allow us to unify the results. Again, we can use a time average 
over iterates of uncoupled Chebyshev maps to generate the weighting with the density po in 
the above integrals. 

3.4 Further iteration steps 

As one can see in Fig. [121 the first iterate of our convolution scheme as given by Eq. ( |T9l) already 
fits the invariant density reasonably well. However, it does not reproduce the structure near 
X = 16 ... 17. Since T2 has an unstable fixed point at = 1 with slope +4, this structure is 
a consequence of local iterates of the CML passing the previously described correlated border 
area one step earlier. 

In a recursive way, we thus define density contributions (iterates of our convolution scheme) 
p^aK^ ~ '^^) (P > 1) follows: 
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• perform the backstep transformation 

• calculate the two-point function p^^~^)(l — ax\ 0) by integrating over 0'_ 

• normalise by p'^^~^\l — ax') 

• integrate over x' (the normalisation cancels) 

• retract the backstep transformation 
Following this recipe, we obtain 

(v)n ^ f Po(0+)^0+Po(0-)^0- rP(A).= lYlkMt:l 

AP7r^2a{l - a) J ^x/Ap + rf (0+) + ' ' 2 ^=o 

(20) 

The entire invariant density is given by the sum over all these contributions, 

p.(i-«x)-Epi"(i-«x) = E /Po(0.M^.M^-W pi) 

p.i ,.i4P,ry'2a(l-a)^ ^x/4>' + (^,+) + r?(0_) 

The series converges rapidly, so usually it is sufficient to take into account the first few functions 
p^^^ only. Roughly speaking, the contribution p^^^ takes care of those points that passed the 
nonhyperbolic region p steps ago. 

Apparently Eq. fl20|) does even hold for the case p = 0, taking into account that the integral 
is invariant under the replacements 0± —<P±- However, in this case the left hand side has to 
be p^^\ax — 1), i.e. the zeroth iterative of our scheme located close to = —1. 

Fig. [12] shows our analytic results ^^^■'(1 — ax), p^a\^ ~ '^^) ^^cl the sum of these two 
contributions. There is excellent agreement with the numerical histogram. The agreement is 
so good that the sum of the two terms is not visible behind the numerical data points. 



4 Explaining the selfsimilar shape 
of the density patterns 

It is obvious that the square root in the denominator of Eq. fl2T]) has essential infiuence on the 
shape of the invariant 1-point density of the CML, since the argument of the square root needs 
to be positive which determines the integration range in our formulas. Abrupt changes of the 
density pattern are connected to zero-border crossings of the arguments of the square roots. 
Looking more closely at the function 

(0+, 0_, x) = ^ + rP(0+) + r^(0_), (22) 

we can find reasons for the patterns. In order to do this analysis, we generated the contour plots 
of r2(0+) + rf (0_) for the values p = 0, 1, 2, 3 by using MATHEMATICA. These plots are shown 
in Fig. [131 It is obvious that the demand r2(0+,0_,a;) > describes a nontrivial limitation of 
the basic integration range 0± G [—1,1]. The equation r2(0+, 4>-,x) = for different values of 
X describes the analytic curves which appear as contour curves in Fig. [131 
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-2-10 1 2-2-10 1 2 

(a) (b) 




-2-10 1 2-2-10 1 2 



(c) (d) 

Figure 13: Contour plots of the function r2(0+) + ^2(0-) for the values p — (a), p = 1 (b), 
p — 2 (c), and p = 3 (d) in the range given by 0± e [—2, 2]. Dark shadings indicate low values, 
bright shadings high values of the function in [—2, 2]. 
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4.1 Patterns generated by p = 

The analytic curve given by the vanishing of the radical 

r°(0+, y)=y + r°(0+) + r°(0-) = y + (0+ + <P-)/2 - 1 (23) 

describes an off-diagonal straight line with intercept 2{y — 1) that limits the integration range 
to lower values of 0+ and (/>_ (cf. Fig. [TST a)). We observe that for |/ = the basic integration 
range is totally suppressed, the integration area then increases like (2j/)^/2 up to y = 1 and like 
(4+ (2 — ?/)^)/2 up to ?/ = 2. Starting from y = 2, the basic integration range is no longer limited 
by the contraint ^2(0+, y) > 0. In accordance to this, the zeroth-order density contribution 
starts with the value zero at y = 0, shows non-differentiable cusps at y = 1 and y = 2 and 
monotonic behaviour elsewhere (compare Fig. [9]). 

4.2 Patterns generated by p = 1 

Here we need to analyse the radical 

ri(0+, 0_, x) = I + rli^+) + r^(0_) = 1 (x + (0+ + 1)^ + (0_ + l)^ _ s) , (24) 

For a; < 8 the analytic curve parametrized by rl{(j)^, x) = is a circle in the (0+, (/)_)-plane 
with radius \/8 — x and center located at (—1,-1) (cf. Fig. [TST b)). The intersection of the 
outside of this circle and the inside of the rectangle — 1 < 0± < 1 determines the integration 
range. 

• For X = the rectangle is totally covered by the circle. In this case, the integration range 
vanishes. We see that the p = 1 density contribution vanishes for x = as well. 

• Between x = and x = 4 the circle uncovers the two positive edges of the rectangle, given 
by the lines 0± = +1. We see that the p = 1 density contribution grows monotonically 
in this interval. 

• Between x = 4 and x = 8 the circle uncovers also the two negative edges of the rectangle, 
given by the lines 0± = — 1. In this interval the p = 1 density contribution declines. 

• Above X = 8 the whole rectangle is uncovered, so there is no restriction of the integration 
range anymore. Starting from x = 8 there are no further cusps. Caused by the increasing 
denominator, the p = 1 density contribution declines. 

4.3 Patterns generated by p = 2 

The second iterate p^\l — ax) shown in Fig. [12] has a cusp at x = 17 which is responsible for 
the shift of the cusp of the invariant density by an amount of — log4(17/16) in Fig. [61 In a 
similar way as for p = 1, we can proceed for p = 2. The only point is that we do no longer 
have a name for the algebraic curve described by r2(0+,0_,x) = 0. However, the behaviour 
is similar as the shape of the curve is a kind of mixture between a circle, a rectangle and a 
triangle (cf. Fig. [TST c)). If we analyse the analytic curve for = 2 more closely, we find that 
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-2-10 1 2 

Figure 14: the function r2(0) for the values p = 0, 1, 2, 3 

the positive edges of the rectangle constituting the basic integration range are fully uncovered 
not at X = 16 as could be assumed for the circle but at x = 17. However, starting from 
X = 16 the positive edges start to be uncovered from the negative side. Therefore, we obtain 
the two cusps at X = 16 and x = 17 as observed. Finally, the basic integration range becomes 
totally uncovered for a value x ~ 33.524, giving rise to a last cusp of this pattern at distance 
— log4(x/16) ~ —0.534 from the first cusp (corresponding to x = 16). 

4.4 Patterns generated by p > 3 

If we look at the function r2(0) for even higher values of p we make an iteresting observation. 
While for the outside region |0| > 1 the function r2(0) starts to diverge, it rapidly converges 
to a limiting function inside the region |0| < 1. This behaviour is shown in Fig. [HI Only the 
inside region is essential for the integration. Therefore, we can conclude that r2(0) ~ 
for relatively small p and |0| < 1. As a consequence of this, there is invariance of the patterns 
generated by the higher iterates of our scheme under the replacement x ^ x/4. As a general 
observation for the radical function r2(0) we add that = +1) = 0, rf(0 = 0) = — 3/4 (for 
p > 0), and 7*2 (—0) = ^2(4') ~ which is caused by the symmetry property of the Chebyshev 
polynomial T2(0) and its iterates T2q{4>) = T2(T29-i (0)). 

5 The scaling behaviour for coupled 
3rd order Chebyshev maps 

All the steps we have performed in order to explain the scaling behaviour of the 2B-string can 
be performed for the 3B-string (coupled 3rd-order Chebyshev maps with s = 1,6 = 0) as well. 
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Scaling regions of the density are transferred by convolution, 



Pa{<P') = - I Pa{<P)Paa {{<P' " (1 " a)T,m / o) d<P. (25) 



However, in order to integrate over the argument (p of the density Paa{4') we have to solve the 
cubic equation 



The three real solutions are 



403 _ 30 = r ^ =: 0. (26) 

1 — a 



cos ( ^ arccos0^ , +^ < 0o < +1 



1 7 27r\ 1,1 

-arccos0+yl, -- < < +- 

1 ^ 47r\ . , 1 



0_ = COS arccos0 + — j , -1 < < -- (27) 

describing the three branches of the inverse map T^^{(j)) (similar formulas apply, of course, to 
the A^-th order Chebyshev map, where there are real solutions). In order to replace the 
differential d(f) in Eq. (125!) by c?0, we have to specify which of the branches are involved. The 
demand |0| < 1 leads to 

0' — a - 0' + a , , 

r < < 11^. 28 

5.1 Iterates close to = +1 

If we use 0' = 1 — ax, is restricted by 

1 — a(l + x) - 1 + a(l — x) , , , , 

1-ax^ ^ < < ^ ^l + a(2-x) (29) 

1 — a 1 — a 

and by— 1<0<+1. turns out to be almost +1, its pre- images under the Chebyshev map 
= T3(0) are found at = +1 and = —1/2. The unperturbed starting contribution of our 
convolution scheme, therefore, comes from = —1/2, involving the branches 0o and 0_. We 
obtain (for 0^1) 

"00 = — sm I - arccos H I — , ~ — =, [Su) 

^3 3 y 3 /^ _ l2 2 v/3(l - 02) 



sm -arccos0+— — — . (31) 

^3 3 y 3 _ l2 2\/3(l - 02) 



Starting with the integration range for as indicated before, both integrations over 0o and 0_ 
are mapped onto this integration range. If we take into account that the integration range for 
00 declines from = — 1/2 to lower values, we can combine both d(j)Q and (i0_ to obtain 

—dcj) ad(f) ad(f) ^^^^ 



3(1-02) ^i^a)^?>{l-(l-ax-a<pY/{l-aY) \jQa{x-l + 

20 



where we inserted ^ for 0' = 1 — ax and calculated the leading order approximation in a. Again, 
the right hand side of Eq. fl25p needs to be calculated only to leading order in a. Therefore, we 
replace all density indices on this side by 0. Finally, as is close to = —1/2, we can replace 
Po(0) by the constant value pq{—1/2) = 2/n^/3 to obtain the zeroth iterate 

p(°)(l-aa;) = ^ r poo{4>)d<f) ^ 2 /• 2po(20 - 0_)po(0_ 
2 f po(0+)'^0+Po(0-)# 



3nV2aJ - 1 + (0+ + 0_)/2 



-) (33) 



The zeroth-order two-point density function p^^^ (0, 0) is obtained by the backstep transforma- 
tion 

The first iterate of our convolution scheme is again the integral of the two-point function, 

p^^\l-ax)= /pW(l-ax',0')rfx' (35) 



where 

0' = - ((1 - ax) - (1 - a)r3(l - ax')) ^ 9x' + 1 - x. (36) 
a 

Integrating over 0' instead of x' and substituting 0' = (0'^ + 0'_)/2, we can again revert the 
backstep transformation to obtain the first iterate 

p«(l - ax) = V / Po(0^)c^0^Po(0-)^0- ^3^^ 



277rV2a J + (0+ + 0_ - 2)/2 + (T3(0+) + T3(0_) - 2)/18 

Continuing in the same manner, we obtain 



00 



P=o p=oQP3nV2aJ Jx/9p + (0+) + r^(0_) 



where 

rU4') = lt^^^. (39) 

2 q=0 y 

5.2 Iterates close to = — 1 

We no longer have to go into details for the calculation close to the lower boundary = — 1. 
The unperturbed domain for the zeroth iterate is given by a region close to = +1/2, the 
higher iterates remain close to = — 1. The final result reads 



00 



^ Z._. ATT / ^,/9. + rf(0+)+rf(0.) 



p=o p=o 9P37rV2a 
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-2-10 1 2-2-10 1 2 

(a) (b) 




-2-10 1 2-2-10 1 2 



(c) (d) 

Figure 15: Contour plots of the function r^{(p+) + rf (0_) for the values p — (a), p = 1 (b), 
p — 2 (c), and p = 3 (d) in the range given by 0± e [—2, 2]. Dark shadings indicate low values, 
bright shadings high values of the function in [—2, 2] 
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-2 

-2-10 1 2 

Figure 16: the function rf (0) for the values p = 0, 1, 2, 3 



We obtain Pa{ay — 1) = Pa(l — clu), i-e. there is symmetry of the 1-point density as expected 
from general principles . Finally, we add the remark that the zeroth iterate for the 3B-string 
close to the lower boundary = — 1 should coincide (up to a general factor) with the zeroth 
iterate (and, therefore, with the density) of the 2B-string close to = — 1, see section 5.4 for a 
numerical check of this fact. 



5.3 The SB-pattern 

What to predict about the 3B-pattern? In Fig. [15] we display the contour plots of the function 
r^{4>+) + r^{4>-) for the values p = 0,1,2,3. There is a step in passing the off-diagonal for 
|0±| > 1) whose size increases with p. However, in the region relevant for the integration 
(— 1 < (f)± < +1) the function converges. In Fig. [16] we show rf(0) for the same values 
p = 0,1, 2, 3. Up to a shift which is given by 

^p(Q\ _ ^ ^ _ ^ ~ (V9)P+^ p^op 9 

^ 2^99 2(1-1/9) 16' ^ ^ 

the odd symmetry of the Chebyshev polynomial T3(0) and its iterates T^g^cj)) = T3(T3q-i(0)) 
cause r^{4>) — r^{0) to have odd symmetry. On the other hand, as in case of the 2B-string we 
have r3(0 = +1) = because of T^q^cp = +1) = +1. Taking these two facts into account, we 
can determine the values of r^{(f)) at = — 1 to be two times the value at = 0. Analysing 
finally the inequality + r^{(j)j^) + r3(0_) > that determines the restriction of the basic 
integration range — 1 < 0± < +1, we can understand the position of the cusps of the density 
curve. 

• For p = the restriction reads x — 1 + (0+ + 0-)/2 > 0. The pattern turns out to 
be the same as for the 2B-string. We obtain full cover for x = (i.e. full suppression 
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of the integration range), half cover for x = 1, and the disappearence of the covering 
triangle from the basic (rectangular) integration range for x = 2. This reproduces the 
three observed segments of the density curve. 

• For p = 1 the restriction reads x — 10 + 2(j)^ + 30+ + 2(f)^ + 3(/)_ > 0. Therefore, we should 
obtain half cover at a; = 10 and disappearence of the cover for x = 20. 

• For high (p-th) iterates, we can apply the limit in Eq. (HTI) . We observe that the half 
cover should be reached for x = (9/8)9^, the disappearence of the cover for x = (9/4)9^. 
All these numbers determine non-differentiable points of the density. 

5.4 Comparison with numerical results 

As we did in Sec. 2, for numerical purposes it is useful to analyse the adjusted energy difference 
a~^/^Po(0)^e(0) as a convenient observable for which we can visualise the scaling behaviour 
in (f)* = — log9(l — (p) close to (p = +1 (because of the symmetry of the spectral density 
Pa(0) in case of the SB-string, the observation of the region close to = — 1 gives the same 
result and, therefore, is obsolete). In Fig. [T7] (corresponding to Fig. IHof the 2B-string) we show 
logarithmic plots of a-^/Vo(0) Ae(0) for the three values a = 0.0000476025, a = 0.00000528917, 
and a = 0.000000587685 of the coupling. The various iterates of our convolution scheme are 
shown in Fig. [18], which corresponds to Fig. [61 We observe 

• for p = two cusps at — logg(ax) = — log9(a) and — log9(2a) 

• for p = 1 the shift of these cusps by an amount of — log9(10/9) = —0.048 . . . 

• for p >2 the shift by an amount of — log9(9/8) = —0.054 . . . 

Finally, in Fig. [19] we compare the numerical result for the invariant density pa(l — ax) (solid 
curve) with the zeroth iterate p^°-*(l — ax) (dashed curve), the first iterate p^^''(l — ax) (dotted 
curve), and the sum of both (dashed-dotted curve). The accuracy of the agreement is again so 
high that one is not able to distinguish between the solid and the dashed-dotted curve. 

6 Proof of the scaling behaviour in a 

At the end of this paper, having all the necessary results at hand, let us finally prove that 
scaling in implies scaling in a. To be more definite, we prove in the following that if V{(f)) is 
a generic test function, then the difference of expectation values taken with invariant densities 
for small coupling a and vanishing coupling a = scales with ^/a and exhibits log-periodic 
oscillations of period logA^^: 

{V{<P))a - (V(0))o = V^fv (log;v^(a)) , fv{x* + 1) = Mx*). (42) 

Examples of this oscillating behaviour were shown in Fig. 1. 
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-log9(a) ^ 



01 2345678 

0* = -|ogg(l -0) 

Figure 17: Adjusted energy difference a~^/^po{^)Ae{(f)) of the SB-string for a = 0.0000476025, 
a = 0.00000528917, and a = 0.000000587685 (from top to bottom). Also indicated is tlie 
breakdown point of the pattern at 4>* = a* = — log9(a) (vertical line). 
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Figure 18: Adjusted energy difference 

a"^^Vo(0)Ae(0) for the SB-string at a = 0.0000476025, 
a = 0.00000528917, and a = 0.000000587685 (from bottom to top). 




Figure 19: 1-point density Pa(l — ax) (solid curve) of the SB-string for a = 0.0000467025 as 
a function of x, compared to the zeroth iterate p^°-'(l — ax) (dashed curve), the first iterate 
p'^^l — ax) (dotted curve), and the sum of both (dashed-dotted curve) 
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6.1 Scaling of pa{l — ax) in a 

To summarise, we found in section 3 that in leading order in a 

P.{l-a.) = -^t I , Ar = 2,3 (43) 

The sum starts at p = for = 3 and at p = 1 for N = 2. Thus we have 

Pa{l - ax) = -/^9{x) (44) 



where the function g{x) is independent of a for small a. 
6.2 Scaling of pa{l — ax) in 

We observed that the adjusted energy difference a~-'^/^po(0)^e(0) scales as 

a-VVo(0)Ae(0) = -^f* (-log^^ (^)) ' ~ = ^^^^ 

close to = +1 (the constant prefactor is chosen for later convenience). If we write 

M,)^.„(,)-e„(,)^^-^. (46) 

we can expand Pa(0) = Po(0) (1 ~ 27r^pQ(</))Ae(0)) to obtain 
Pa{<P) ^ po(0)(l + vrVo(0)Ae(0)' 



- Po(0) + ^P^(0)r(-log^.(^)). (47) 

In order to prove the periodicity property f*{x* — 1) = f*{x*) we use the fact that we can 
represent the invariant density Po{l — ax) by a series expansion analogue to Eq. P3|) . 

, 1 2^1 2^1 



Try^l - (1 - ax)2 N ^f^^, A^p A^ttv^ ptS^i N^p^x/N^p 
2 ^ r po{(j)+)d(f)+po{(f)^)d(j)_ 

' mp^Jx/mp 



NnV2a p=o,i 



In terms of /(x) = /*(— log^2(x)) = f*{x*) the periodicity reads f{N'^x) = f{x). The function 
/(x) is given by 



f^^^ = ^ y f ( Po(0+)ci0+po(0-)c?0_ _ po(0+)t/0+po(0_)rf0, \ = V / 



(49) 
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The periodicity can be shown for values x ^ A^^^ where P again is chosen large enough to 
guarantee the p-independence of the function r^{(f)), r^{4>±) ~ fN{4'±) fo^' P ^ P- For this 
range of values for x we have 



fp{N^x) = N^x 




mP^mx/N^P + r^(0+) + r^(0-) N'^P^N^x/N^P 
Po (0+ ) d(l)+pQ (0_ ) po (0+) rf0+Po (0- ) dcp- 




On the other hand, for p < P the value of x/N'^'^ is large. Therefore, we can expand the 
radicals, 

1 11/ A^^P \ 



2a; 



1 /A^2p\3/2 

= -2(^j (^^('^+)+^^('^-))- (51) 

Because / po(0+)(i0+po(0-)(^0-(r5^(0+) + r^(0-)) remains finite, we can neglect the above 
difference in fp{x) for small values of p. Hence, we have finally shown that 

oo oo oo 

f{N'x) = Y: fp{N'x) ^ E ^ ^ f,ix) = fix). (52) 

p=0,l P=0,1 p=o,i 

We have shown that this equation holds for values x with ax <^ 1. Note that even though it is 
obvious from Figs. H] and [T7| that at about ax ^ 1 the pattern is still periodic, we cannot prove 
this by using the analytical formula fH51) . which is a perturbative result valid close to the upper 
border of the interval. Also, the approximation we did in Eq. (HH!) no longer holds. However, 
if a is small enough, then the region in question constitutes only a small fraction of the entire 
repetitive pattern. 

6.3 Scaling of expectation values in a 

Let us now prove the scaling behaviour of 

(V^(0)), - {Vmo = £^ {Pa{<f>) - Po(0)) ^(0)#, (53) 

where ^ is a generic test function of the local iterates of the CML. In order to apply the two 
scaling relations of the previous subsections, we split the integral over into two parts by 
chosing a subdivision point = 1 — ca. c is arbitrary but fixed (for N = 2 e.g. c = 2^ = 128). 
We obtain (V^(0)), - (V^(0))o = 4(c) + hie) where 

Iu{c) = C (p,(0)-po(0))V(0)rf0, (54) 

Jl—ca 

/I— ca 
(Pa(0)-Po(0))V^(0)#. (55) 
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The fact that the integral /„(c) scales with ^/a can be shown very easily. From Eq. (jH]) we 
obtain 



c 



/u(c) = a (pa{l — ax) — pQ(l — ax))V(l—ax)dx 
Jo 

a r l^g(x) ^ I V(l - ax)dx ^ ^/aV(l) T I g(x) ^ 1 dx (56) 

Jo \va Tw/lax) Jo \ 



Try 2a; ^ 

where the approximations hold for ax <^ 1. If we replace a by Aa where A is (so far) an arbitrary 
factor and c is kept constant, we obtain 

4(c) ^ v^r(l) r L(x) ^ ) dx = VXIuic). (57) 

"'0 \ 7rv2x / 

For the second integral Ii{c) we can assume that is large enough to guarantee the preser- 
vation of the pattern given by Eq. (1521) . We use Eq. (H7I) to obtain 

/,(c) = v^/7r(-.og,„(i^))nw, y(4>)-j^^^y (58) 

In replacing a by Aa we obtain 

rl—\ac I / ][ 



li{c) ^v^J^r log^2 [^-^j + log7v^(A)j V^(0)#. (59) 

In general, we do not have /;(c) \/~XIi{c) as in the case of the integral Iu{c). However, if we 
chose A = we can use that f*{x*) is (nearly) periodic with period 1 to obtain 

h{c) - -^a /* (^-log^. ( a j j ^^'^^^'^ = ]^^'(c/iV')- (60) 

If we replace the upper limit of the integral by 1 — ca, we neglect an integration range in of 
length of order 0(a). Therefore, in leading order we obtain J;(c) Ii{c)/N. Collecting both 
results, we finally obtain that under the transformation a — > a/N"^ we have 

{vm, - (v(0))o - {vm,,^. - {v{4>))o = ^((^(0))a - {vmo) m 

whereas for general transformations a Aa we obtain a result different from a simple multi- 
plication by a factor V^- It is obvious, then, that we have 

{Vm, - (V(0)o = v^/;; (log;v^(a)) (62) 

where /]^(log^2(a)) is a periodic function with period 1 depending on the test function V^(0). 

Conclusion 

The treatment of nonhyperbolic CMLs exhibing chaotic behaviour is notoriously more difficult 
than that of hyperbolic ones. In this paper we tackled the problem by introducing an iterative 
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convolution technique that provides a perturbative expression for invariant 1-point and 2-point 
densities. For small coupling parameters a the invariant 1-point density Pa of the CML is 
obtained as a sum of density contributions where the index p essentially decribes at which 
time in the past the nonhyperbolic region is active in our convolution scheme. We obtained 
explicit perturbative expressions for the densities of diffusively coupled Chebyshev maps of N- 
th order. Our main examples were the cases N = 2 and N — 3, which significantly differ in 
their symmetry properties. While in the uncoupled case the invariant 1-point density is smooth, 
an arbitrarily small coupling a induces a selfsimilar cascade of patterns with a variety of cusps 
and non-differentiable points, which can all be understood by our perturbative approach. We 
proved that arbitrary expectation values scale with the square root of the coupling parameter 
and that there are log-periodic oscillations both in the phase space and in the parameter space. 
Though, for reasons of concreteness, we mainly dealt with 2nd and 3rd order Chebyshev maps, 
our techniques can be applied in a similar way to other nonhyperbolic coupled systems where 
the local maps exhibit fully developed chaos. 
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